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We present the results of our theoretical study on the effects of geometrical frustration and 
the interplay between spin and orbital degrees of freedom in vanadium spinel oxides AV2O4 
(A = Zn, Mg or Cd). Introducing an effective spin-orbital-lattice coupled model in the strong 
correlation limit and performing Monte Carlo simulation for the model, we propose a reduced 
spin Hamiltonian in the orbital ordered phase to capture the stabilization mechanism of the an- 
tiferromagnetic order. Orbital order drastically reduces spin frustration by introducing spatial 
anisotropy in the spin exchange interactions, and the reduced spin model can be regarded as 
weakly-coupled one-dimensional antiferromagnetic chains. The critical exponent estimated by 
finite-size scaling analysis shows that the magnetic transition belongs to the three-dimensional 
Heisenberg universality class. Frustration remaining in the mean-field level is reduced by ther- 
mal fluctuations to stabilize a collinear ordering. 
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1. Introduction 

Spinels are one of the most typical geometrically- 
frustrated systems. Among them, the so-called B spinel 
oxides AB2O4, where A cations are nonmagnetic, have 
attracted much interests since magnetic B cations 
form the pyrochlore lattice, which consists of a three- 
dimensional (3D) network of corner-sharing tetrahedra. 
The system suffers from strong magnetic frustration: 
Magnetic correlations are strongly suppressed and com- 
plex phenomena may appear at low temperatures due to 
a large number of nearly-degenerate ground states. 

In this paper, we will investigate typical B spinels 
AV2O4 with divalent nonmagnetic A cations such as 
Zn, Mg or Cd. Besides the geometrical frustration, these 
vanadium spinels have another key issue, the orbital de- 
gree of freedom. Since each V 3+ cation has two 3d elec- 
trons in threefold ti g levels, the system has the orbital 
degree of freedom in addition to the spin. Hence, vana- 
dium spinel oxides are intriguing systems which provide 
two important issues in strongly correlated systems, i.e., 
geometrical frustration and the interplay between spin 
and orbital degrees of freedom. 

Vanadium spinels AV2O4 (A — Zn, Mg or Cd) show 
two different phase transitions at low temperatures; one 
occurs at around 50K and the other is at around 40K. 1 
Note that the transition temperatures are significantly 
lower than Curie- Weiss temperature ~ 1000K, 2 which is 
attributed to the effect of geometrical frustration. The 
former is a structural transition from high-temperature 
cubic phase to low-temperature tetragonal phase. The 
latter is a magnetic transition with a collinear antiferro- 
magnetic (AFM) ordering. The magnetic structure con- 
sists of the staggered ordering with T-|-T _ |-' ■ ■ structure 
along the xy chains and the period-four ordering with 
T-T-|-|-" ' ' structure along the yz and zx chains 3 [see 
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Fig. 1 (b)]. (We take the z axis in the tetragonal c direc- 
tion.) The issue is the microscopic mechanism of the two 
transitions: How is the degeneracy due to the geometrical 
frustration lifted? What is the role of the orbital degree 
of freedom? 

The authors have proposed an effective spin-orbital- 
lattice model for this problem and shown that the two 
transitions are well reproduced by Monte Carlo (MC) 
simulation performed for this model. 5 ' 6 In the present 
study, we will focus on the magnetic frustration in the 
spin-orbital coupled system. We will introduce a reduced 
spin model to explain the low energy physics in the or- 
bital ordered phase, and clarify the nature of the mag- 
netic transition including its critical properties. 

This paper is organized as follows. In Sec. 2, we intro- 
duce the effective spin-orbital-lattice model which has 
been derived in the previous publications, 5 ' 6 and de- 
scribe the method of calculations. In Sec. 3, we propose 
a reduced spin Hamiltonian in the orbital ordered phase 
and explain the microscopic mechanism of stabilizing the 
complex AFM ordering. Numerical analysis of the critical 
exponent is also shown. Section 4 is devoted to summary. 

2. Model and Method 

We start from a multiorbital Hubbard model with 
three t2 g orbitals on the pyrochlore lattice, and consider 
the perturbation in the strong correlation limit to de- 
scribe the low energy physics of the insulating vanadium 
spinels. By using atomic eigenstates with two 3d elec- 
trons per site in a high-spin state as the unperturbed 
states, we derive the effective spin-orbital-lattice coupled 
model (the so-called Kugel-Khomskii type Hamiltonian 4 ) 
in the form 5 ' 6 



H — H. n 



H 



JT ■ 



H so — — J hij — J3 hij, 

m «<j» 



(1) 

(2) 
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hij — {A + BSi • S j)[ <l ^ia(ij) n ja(ij) n ia(ij) ^ja(ij) ] 
+ C(1 — Si ■ S j)^i a (ij) ra ja(ij); (3) 

ffjT = 7E^ e * + E^/ 2 - A E^' ( 4 ) 

i i (ij) 

where Si is the S — 1 spin operator and ni a is the den- 
sity operator for site i and orbital a = 1 (d yz ), 2 (d zx ), 
3 {d xy ). Here, ¥ii a = 1 — ni a and we impose a local con- 
straint ^ a=1 n^Q, = 2 at each site. The summations with 
(ij) and ((ij)) are taken over the nearest-neighbor (NN) 
sites and third-neighbor sites, respectively. Here, we take 
into account only the dominant cr-bond hopping integrals 
in the original multiorbital Hubbard model, for instance, 
hoppings between d xy orbitals in the same xy plane. 5,6 
This approximation results in the orbital diagonal inter- 
action in H so ; a(ij) is the orbital which gives rise to the 
a bond between sites i and j. Hjt describes the orbital- 
lattice coupling part, where 7 is the electron-phonon cou- 
pling constant of the tetragonal Jahn- Teller (JT) mode, 
Qi denotes the amplitude of local lattice distortion at site 
i, and e; = riji+n^ — 2^3. A describes the interaction be- 
tween NN JT distortions, which mimics the cooperative 
aspect of the JT distortion. 

An important feature of the model (1) is the symmetry 
of the intersite interactions. In the spin part, the inter- 
action is isotropic Heisenberg type. On the other hand, 
in the orbital part, it depends only on the density op- 
erator, that is, there is no transverse component which 
mixes different orbitals. Hence, the orbital interaction is 
three-state Potts type corresponding to three t2 S states. 
Moreover, it depends on bond direction and also orbital 
states at both ends of a bond. These peculiar proper- 
ties of the orbital interactions are crucial to trigger the 
orbital ordering by lifting the degeneracy in the present 
frustrated system. 5 ' 6 

The parameters in H so are given by the coupling con- 
stants in the starting multiorbital Hubbard Hamiltonian, 
and in the following we use the reasonable estimates 
given in Ref. 6 as J 3 /J = 0.02 with J ~ 200K, A = 1.21, 
B = 0.105, and C = 0.931. For the JT parameters, we 
take 7 2 /J = 0.04 and A/ J = 0.15, which are typical val- 
ues to have the tetragonal distortion consistent with the 
experimental result. 6 Hereafter, we will set the lattice 
constant of the cubic unit cell as a length unit and use 
the convention of the Boltzmann constant fcg = 1. 

We have studied thermodynamic properties of the 
model (1) by employing classical Monte Carlo (MC) sim- 
ulation to avoid the negative sign problem due to the 
geometrical frustration of the pyrochlore lattice. Since 
quantum nature exists only in the spin 5 = 1 operators 
in H SOl we approximate them by classical vectors with 
unit modulus. 6 We typically perform 10 5 MC samplings 
for measurements after 10 5 steps for thermalization. Sys- 
tem sizes are up to L 3 = 12 3 in the cubic units, which 
includes 12 3 x 16 vanadium sites. 

3. Results 

3.1 Reduced spin model under orbital ordering 

MC results show that as temperature decreases, first 
an orbital ordering takes place with the tetragonal JT 




Fig. 1. (a) Effective spin exchange interactions for nearest neigh- 
bors in the orbital ordered phase in a cubic unit cell of the py- 
rochlore lattice. White (gray) circles are the V sites where the 
d yz and d xy (d zx and d xy ) orbitals are occupied. Solid (dashed) 
bonds represent the strong AFM interactions Jaf (the weak FM 
interactions Jp) in the xy (yz and zx) planes, (b) Effective spin 
exchange interactions for third neighbors, Jg. Two sublattices 
are shown which are connected by Jg. 

distortion, corresponding to the structural transition in 
experiments. 6 This is triggered by the highly anisotropic 
nature of the three-state Potts-type orbital interaction 
described above. This transition is discontinuous, and 
the orbital moment shows a large jump at the transition 
temperature To — 0.19J. The ordering structure is A 
type, i.e., the layered order with alternative stacking of 
(dzx, d xy ) and (d yz , d xy ) occupied planes in the z direc- 
tion as shown in Fig. 1 (a). 

Since the orbital polarization below To is large and 
quickly converges to the saturated value, 6 it is conve- 
nient to consider a reduced spin Hamiltonian by freezing 
the orbital degree of freedom in order to capture the low 
energy physics in the orbital ordered phase. The reduced 
spin Hamiltonian is obtained by replacing the orbital 
parts in Eq. (2) by their mean values in the fully po- 
larized state. For instance, we replace nan^ and nafij^ 
in H so by (n i3 n j3 ) = 1 and (n i3 h j3 ) = 0, respectively, 
for the bonds in the xy planes. Then we end up with the 
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reduced model in the form 



, J/^pSi-Sj- 

(ij)exy 



{ij)£yz,zx {{ij)) 

(5) 

The first two terms describe the NN spin exchange in- 
teractions, in which the former (latter) summation is 
taken over the NN pairs in the xy (yz and zx) planes. 
Jaf = JC > is antiferromagnetic (AFM) and Jp = 
— JB < is ferromagnetic (FM). These exchange inter- 
actions are shown in Fig. 1 (a). The last term describes 
the third-neighbor interactions, where the summation is 
taken over the third-neighbor er bonds shown in Fig. 1 
(b). J' z = J 3 C > is AFM. Here we omit the orbital- 
lattice part for simplicity. 

First we consider only the NN exchange interactions. 
One important point is the magnitude of two different in- 
teractions Jaf and Jp. By using the parameters in Sec. 2, 
we obtain Jaf — 0.931J and \Jp\ — 0.105J, and hence 
the AFM exchange in the xy planes is about ten times 
stronger than the FM exchange in the yz and zx planes. 
Thus, the system can be regarded as one-dimensional 
(ID) AFM chains coupled by weak FM exchange inter- 
actions. Another key feature is the frustration due to the 
interchain coupling Jp. Once the AFM correlation de- 
velops along the xy chains by the strong Jaf, coupling 
of the xy chains are frustrated; in the mean-field level, 
the total energy does not depend on the relative angle 
of AFM moments in any two xy chains. Therefore, the 
system is reduced to independent ID AFM chains and it 
is hard to establish a three-dimensional (3D) AFM order 
at this stage. 

The frustration remaining among the xy chains is 
almost reduced by the small third-neighbor exchange 
J3 ~ 0.0186J. As shown in Fig. 1 (b), J3 connects the 
parallel xy chains and stabilizes a 3D AFM ordering. 
Note that however there still remains frustration between 
two sublattices; one consists of (d zx , d xy ) occupied sites 
and the other consists of {d yz , d xy ) occupied sites. The 
relative angle between two sublattice moments M\ and 
M2 is free in the mean-field level. Reduction of the frus- 
tration will be discussed in Sec. 3.3. 

The picture based on the reduced spin model (5) is con- 
firmed by our MC calculations. Figure 2 shows the tem- 
perature dependences of the NN and third-neighbor spin 
correlations in the effective spin-orbital-lattice model (1): 

S NN(3rd) = T,(ij)(({ij)))e v (Si ■ Sj)/N h , where N h is the 
number of bonds in the summation and v — xy,yz, zx. 
Note that S^ yz ^ equals to S^ zx ^ by symmetry. Below the 
orbital ordering temperature To, AFM correlations de- 
velop along the xy chains (S^^ < and S^' > 0), 
while correlations in the yz and zx chains remain small. 
Particularly, the NN correlations become almost 

zero because of the frustration in the interchain coupling 
Jp. Below Tn — 0.12J (which will be assigned to the 3D 
AFM transition temperature in Sec. 3.2), rapidly 
grows due to J3. Thus the system becomes highly one- 
dimcnsional below To, and 3D interchain correlations de- 
velop below Tn- 

Consequently, the orbital ordering introduces the spa- 
tial anisotropy in the intersite spin exchange interactions, 
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Fig. 2. Temperature dependences of spin correlations for the sys- 
tem size L = 12. Error bars are smaller than the symbol sizes. 
See the text for details. 
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Fig. 3. (a) System size extrapolation of the staggered moment 
along the xy chains. Lines are the linear fits for the data. Inset: 
System size dependence of Mg. (b) Temperature dependence of 
the extrapolated data in (a). The curve shows the fit by Mg 1 "' oc 
(Tn — T)@ . Error bars are smaller than the symbol sizes. 



which plays a key role to reduce the geometrical frustra- 
tion. The system is effectively reduced to weakly-coupled 
ID chains and finds a way to establish a 3D AFM order. 

3.2 Magnetic transition and universality class 

The 3D AFM order is indeed found to occur in our 
MC results for the model (1). The staggered magnetic 
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Fig. 4. (a) Temperature dependence of the collinearity. (b) Sys- 
tem size dependence. The lines are guides for the eyes. 

moment Ms, which is the magnitude of the staggered 
moment Mi or M 2 , grows continuously atT~0.12J< 
Tq. To determine the transition temperature Tn pre- 
cisely, we perform a finite-size scaling analysis. As shown 
in the inset of Fig. 3 (a), Mg shows non-monotonic 
size dependence within the range of the system sizes 
calculated here, probably because of the large spatial 
anisotropy of the AFM correlations discussed in Sec. 3.1. 
Instead, we use the staggered moment along the AFM 
xy chains, Ma Xy ^ , which is defined by the average of the 
magnitude of AFM moments in the xy chains. 6 The re- 
stricted summation reduces effects of anisotropic corre- 
lations. 

Fig. 3 (a) shows the system-size dependences of Mg . 
The MC data are fitted by a linear function of 1/VZ 
(Ref. 6). The values extrapolated to L — ► 00 are plotted 
as a function of temperature in Fig. 3 (b). 

The transition temperature Tn and the critical expo- 
nent are obtained by the scaling fit M^ v ^ oc (Tn — T) 13 
for the data in Fig. 3 (b). The best fit gives the estimates 
T N /J = 0.1 18 ±0.001 and /3 = 0.36 ±0.01. The exponent 
(3 is consistent with that of the 3D isotropic Heisenberg 
model with short-range interactions, (3 = 0.365 (Ref. 7), 
which indicates that the AFM transition in the present 
system belongs to the 3D Heisenberg universality class. 

3.3 Collinear ordering due to 'order by disorder' 

As explained in Sec. 3.1, in the mean-field level, frus- 
tration remains between two sublattice moments Mi 
and M 2 shown in Fig. 1 (b) even when 3D long-range 
order is well established in each sublattice. Note that ei- 
ther Mg or M { b xy) in Sec. 3.2 does not give information 
about the relative angle between Mi and M 2 . 

Here, we measure the collinearity by Ci 2 = (cos 2 #12), 



where 6*12 describes the angle between Mi and M 2 . Fig- 
ure 4 shows our MC results of C12 for the model (1). C\ 2 
rapidly increases below Tn as shown in the inset. System- 
size dependences of Ci 2 show that the data approach 1 
in the thermodynamic limit below Tn, which indicates 
that the magnetic order is collinear. 

The collinear ordering is understood by the so-called 
order-by-disorder mechanism. 8 Our MC results show 
that the frustration remaining in the mean-field level is 
reduced by thermal fluctuations and the collinear AFM 
state is stabilized. In Ref. 5, the authors have shown 
that quantum fluctuations also favor the collinear state 
to minimize the zero point energy of magnons. It is well 
known that both thermal and quantum fluctuations favor 
a collinear state in many frustrated spin systems. This is 
also the case in the present system. 

The magnetic structure obtained by MC calculations 
is shown in Fig. 1 (b). The ordering pattern is consistent 
with the experimental result by the neutron scattering. 3 

4. Summary 

We have found that the interplay between spin and 
orbital is crucial in the geometrically-frustrated vana- 
dium spinels. The orbital ordering introduces spatial 
anisotropy in the effective spin exchange interactions 
and drastically reduces the magnetic frustration. In the 
orbital ordered phase, the system can be regarded as 
weakly coupled ID AFM chains. By applying the finite- 
size scaling analysis to numerical data, we have shown 
that the AFM transition belongs to the universality 
class of the 3D unfrustrated Heisenberg model. We have 
also pointed out the importance of thermal fluctuations 
to stabilize the collinear AFM state by the order-by- 
disorder mechanism. 
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